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The High Altitude Water Cherenkov (HAWC) Observatory monitors the gamma-ray sky in the 
100 GeV to 100 TeV energy range with >95% uptime and unprecedented sensitivity for a sur¬ 
vey instrument. The HAWC Collaboration has implemented an online flare monitor that detects 
episodes of rapid flaring activity from extragalactic very high energy (VHE) sources in the dec¬ 
lination band from -26 to 64 degrees. This allows timely alerts to be sent to multi wavelength 
instruments without human intervention. The preliminary configuration of the online flare moni¬ 
tor achieves sensitivity to flares of at least 1 hour duration that attain an average flux of 10 times 
that of the Crab Nebula. While flares of this magnitude are not common, several flares reaching 
the level of 10 Crab have been observed in the VHE band within the past decade. With its survey 
capabilities and high duty cycle, HAWC will expand the observational data set on these particu¬ 
larly extreme flares. We characterize the sensitivity of the online flare monitor and outline plans 
for its upcoming deployment. 
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1. Blazar Flares at Very High Energies 

Blazars comprise a class of active galactic nuclei (AGNs) with relativistic jets oriented along 
the line of sight to Earth [1]. More than 50 objects in this class have been detected in the very high 
energy (VHE; E > 100 GeV) band^ Typical blazar spectra exhibit a characteristic “two-bump” 
feature, with the low energy bump, spanning radio to hard x-rays, being interpreted as synchrotron 
emission from a population of relativistic electrons in the jet. The high energy bump, spanning MeV 
energies and above, is believed to originate either due to a leptonic scenario—in which low energy 
photons are inverse Compton scattered by the leptons responsible for the synchrotron emission—or 
due to a hadronic scenario—in which proton synchrotron emission, pion decay, or other hadronic 
interactions dominate the high energy emission (see for example [2] and references therein). 

Blazars are known to be highly variable sources in the VHE energy band, with fluxes varying 
on time scales ranging from minutes to months and increasing by more than an order of magni¬ 
tude over their baseline quiescent states [3]. The study of these flaring states can offer insight into 
the mechanisms responsible for powering the relativistic jets, including the possibility to discrim¬ 
inate between leptonic and hadronic origins of the VHE emission. Eurthermore, the increase in 
high energy event counts during flares provides an opportunity to measure the extragalactic back¬ 
ground light via the attenuation of primary gamma rays as they propagate through the intergalactic 
medium [4], and the timing characteristics of the flares may be used to measure Eorentz invariance 
violation [5] or detect the presence of an intergalactic magnetic field [6, 7]. Identifying flaring 
blazars and obtaining a data set of multiwavelength observations during the flares is therefore an 
ongoing endeavor of great scientific interest. 

Until recently, VHE observations of blazar flares on hour time scales have been conducted 
exclusively by imaging atmospheric Cherenkov telescopes (lACTs), which observe the Cherenkov 
radiation from secondary particles in air showers initiated by high energy gamma rays [8]. The 
current generation of lACTs has collected numerous observations of flaring blazars. Over the past 
decade, a few particularly extreme flares, reaching more than 10 times the flux from the Crab Neb¬ 
ula (the strongest steady source of VHE gamma rays), have been observed [9, 10]. However, since 
lACTs are pointed instruments with limited sky coverage and are likely to continue observing a 
source found to be in a flaring state, the catalog of VHE flares is biased, rendering the determina¬ 
tion of the rate at which blazars flare difficult. 

In contrast to lACTs, VHE survey instruments observe a wide field of view that allows them to 
monitor a large number of sources simultaneously. The High Altitude Water Cherenkov (HAWC) 
Observatory, located at 19° N latitude at an elevation of 4100 m in Sierra Negra, Mexico, is the most 
sensitive VHE survey instrument ever constructed. The detector comprises 300 optically isolated 
water Cherenkov detectors, each of which contains 4 upward-facing photomultiplier tubes (PMTs) 
monitoring 200,000 liters of purified water for the Cherenkov radiation produced by secondary 
particles propagating through the detector volume. HAWC, described in detail in [11], features a 
field of view of 2 sr, an uptime of >95%, and a substantial improvement in sensitivity over its 
predecessors [12, 13], enabling it to conduct an unbiased survey of flaring activity in blazars. This 
capability will greatly enhance the existing data set of extreme flares, providing new opportunities 

*See http://tevcat.uchicago.edu for an up to date list. 
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to study the blazar central engine, as well as to investigate the propagation of VHE gamma rays 
over cosmic distances. 

Due to its location, HAWC achieves maximum sensitivity to sources with declinations close 
to 19°. For sources culminating at an elevation angle of 45° (declinations ranging from —26° to 
64° for HAWC), the sensitivity is reduced by approximately a factor of 10 [14] and falls rapidly 
outside of this range. For the present investigation, we restrict our attention to sources within this 
declination range. We also focus this study on the sensitivity of HAWC to detect flares on hour time 
scales, with the intention to issue alerts as quickly as possible after such flaring states are identified. 
The sensitivity of HAWC to flares on time scales of 1 day or more and initial observations of such 
flares are presented elsewhere in these proceedings [15]. 

2. The HAWC Online Flare Monitor 

The online flare monitor tracks all blazars known to be VHE emitters, as well as a selection of 
nearby sources (z < 1) from the IFHF catalog produced by the Fermi Collaboration [16]. Events 
triggering the HAWC detector are reconstructed in real time by the HAWC online system, after 
which they are processed by the online flare monitor. The events are divided into analysis bins 
defined in terms of Ufrac, the fractional number of PMTs that participate in the event. The parameter 
Ufrac serves as a crude proxy for the total energy deposited by the air shower at ground level. An 
improved estimator for the energy of the primary particle will be used to construct the analysis bins 
in the near' future. 

Data reduction is achieved via a set of selection cuts that maximize sensitivity to the Crab 
Nebula. The selection cuts include the distribution of charge recorded by the PMTs and the qual¬ 
ity of the fits that determine the shower core and angular position on the sky. Separate optimal 
cuts are determined for each analysis bin. Further information about the cuts is available in these 
proceedings [17, 18]. 

The determination of the analysis bins for the online flare monitor is different from that for 
steady sources as presented in [18]. Because we seek transient signals on hour time scales, the 
event rate in each analysis bin must be much larger than 1 per hour. We therefore define analysis 
bins to achieve roughly equal count rates in each bin, and we optimize the cuts to maximize the 
sensitivity to the Crab Nebula. We have tested the online flare monitor with 1, 2, 3,4, and 8 analysis 
bins, and we find that the sensitivity is maximized for the 3-bin case. 

3. Flare Identification 

We define flares as periods during which the flux from a monitored source exhibits a significant 
increase over its quiescent level. Because the detector response depends on the zenith angle (and 
therefore changes as the source moves thr'ough HAWC’s field of view), simply monitoring the 
rate of events from the source is insufficient to identify flares. Instead, we search for changes in 
the signal to background ratio R = Von/otVoff, where following [19] we define Von and Vou as 
the on-source and off-source counts and a as the ratio of on-source to off-source exposure. The 
off-source counts Vou are derived via the direct integration method [20], in which the product of 
the detector rate and the normalized counts in an angular bin of solid angle D/) in local dectector 
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coordinates is integrated over a pre-defined time period, here chosen to be 2 hours. The value of 
Non is determined from the counts over the same time interval and in an angular bin, centered on 
the position of the source, with the same solid angle Q.jj, but in equatorial coordinates instead of 
local detector coordinates. 

Our search employs the Bayesian block algorithm [21], which we use to identify the optimum 
partitioning of the data into blocks within which R is constant. The edges of the blocks are then 
taken as triggers indicating the start or end of flaring activity. For a given set of Np data points, 
the algorithm finds the optimum partition in Np time, compared to the brute force method, which 
scales as 2^'’. Optimality is defined by a fitness function F that specifies the fitness of a given block 
to be represented by a constant value. The fitness function must depend only on the data within the 
block, and it must be additive such that the fitness of a set of blocks is the sum of the individual 
blocks’ fitnesses. 

The HAWC events are too numerous to permit an unbinned analysis, so we bin the observa¬ 
tions into data points of 1 minute duration. Each data point i consists of a number of on counts, 
denoted by and a number of off counts, denoted by m,. Since the detector response changes 
with time as the source elevation changes, we assume that each 1-minute interval has a separate 
true Poisson mean characterizing the distribution from which the off count rate m, is derived. We 
denote these mean values by /i,. Our hypothesis for each block is that the ratio R is constant, so we 
can write the true Poisson mean of the distribution from which the on counts are sampled as a/r,7?. 
We then follow Section 3 of [21] and construct the fitness function for the block as the sum of the 
log likelihoods for each data point i in the block: 

lnL(H,m|/?,/I) = ^[—a/r,7?-|-n;ln(a/r,7?) —/x, -I- m,Tn/r,-— g(n,■,»?,■)], (3.1) 

i 

where g{ni,mi) = lnr(n,' -|- 1) -|-lnr(m,- -|- 1) and r(v-|- 1) = x\ is the factorial function. To derive 
the fitness function, we choose the values of R and /i,- that maximize the likelihood, obtaining 

+ [('^! + "h) ln(n,- -h nii) - m,-)]. (3.2) 

i 

Here, N = and M = and all sums run over the data points in the block. Since the HAWC 
data are divided into analysis bins as well as time bins, the index b in Equation 3.2 indicates that 
this is the fitness contribution from one analysis bin. The total block fitness is then given by 


Fb = -N 
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F = l,Fb. 

b 


(3.3) 


The final input to the Bayesian block algorithm is a probability p that accounts for our prior 
belief in the number of blocks. Specifically, we assume that it is p times less likely for the data 
to be described hy Nb + I blocks than by Nb blocks, for Ag G {1,2,3,...}. The value of p is then 
tuned to set a predetermined false positive rate that depends on the monitored source class. Eor 
this study, we fix the prior to obtain a false positive rate of 10^^ events per minute per monitored 
source. This requires different priors to be set for different numbers of analysis bins. Eor example, 
the prior in the case of 1 analysis bin is p\ « 10^^, while the prior for 8 analysis bins is pg « 10^^. 
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Figure 1: (a): Example of a simulated flare with a flux 10 times that from the Crab Nebula, with 
the rising edge identified by the flare monitor. The flare monitor was run with only 2 analysis bins 
to limit the number of plots, (b): Probability for the optimal 3-bin flare monitor to detect flares of 
varying strength and duration, as a function of buffer size. We consider only those flares for which 
some part of the flare occurs when the event rate is at least half of its value at source culmination. 
The plot is made for a false positive rate of 10^^ triggers per minute per monitored source. 

4. Flare Monitor Sensitivity 

In order to characterize the sensitivity of the flare monitor, we inject into the observed data a 
number of simulated flares with spectra scaled to the spectrum of the Crab Nebula as measured by 
HAWC. We then process the simulated results using the Bayesian block algorithm from Section 3, 
determining the number of simulated flares that are detected. A sample simulated flare with a flux 
of 10 times that of the Crab Nebula appears in Figure la. For clarity, the figure depicts the results 
from running the flare monitor with 2 analysis bins instead of the optimal value of 3. The flare 
monitor identifies the beginning of the flare as a significant change point, demonstrating its ability 
to detect flares. 

Since the goal of the online flare monitor is to issue rapid alerts, we do not need to run the 
algorithm over the entire cumulative data set. Instead, we run it over a buffer of configurable 
duration that contains only the most recent events. Figure lb depicts the probability for the flare 
monitor to identify flares of varying strength and duration, as a function of buffer size. In this case, 
we have run the flare monitor in the optimal case with 3 analysis bins. As shown in the figure, in all 
cases the detection rate rises quickly at small buffer sizes, but then begins to level off, suggesting 
that buffer sizes larger than 10 hours have little effect on the sensitivity. For a sufficiently large 
buffer, the probability to detect a flare that reaches 10 times the flux from the Crab Nebula for a 
period of 1 hour is approximately 80%. The figure also shows that the detection probability depends 
more strongly on the strength of the flare than its duration: a factor of 3 reduction in duration results 
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in a detection probability of approximately 35%, whereas dropping the flare strength by only a 
factor of 2 produces a detection probability of less than 20%. 

5. Alert Expectations and Future Improvements 

At present, the HAWC online flare monitor alerts are being tested internally by the HAWC 
Collaboration. When the design of the system is finalized, alerts will be made available via email 
when significant flares are detected. We expect to begin issuing these alerts soon. Additionally, 
although the existing monitor already achieves sensitivity to flares that are known to occur, several 
improvements to the HAWC sensitivity are currently being investigated. Among these are better 
gamma-hadron separation techniques and the development of a refined energy estimator to define 
the analysis bins. When complete, these improvements will allow us to issue alerts for flares of 
lower strengths, as well as to more rapidly identify periods of extreme flaring activity. Finally, 
we are developing two additional versions of the online flare monitor, one to search for Galactic 
transients, and another to conduct an all-sky search for transient VHE phenomena. The main 
difference between the blazar flare monitor and these additional monitors will lie in setting the 
Bayesian prior, as rare events or large numbers of sources require a much lower false positive rate. 
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